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Spatial Analysis of Opportunistic Downlink 
Relaying in a Two-Hop Cellular System 
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Abstract 

We consider a two-hop cellular system in which the mobile nodes help the base station by relaying 
information to the dead spots. While two-hop cellular schemes have been analyzed previously, the 
distribution of the node locations has not been explicitly taken into account. In this paper, we model 
the node locations of the base stations and the mobile stations as a point process on the plane and then 
analyze the performance of two different two-hop schemes in the downlink. In one scheme the node 
nearest to the destination that has decoded information from the base station in the first hop is used as 
the relay. In the second scheme the node with the best channel to the relay that received information 
in the first hop acts as a relay. In both these schemes we obtain the success probability of the two hop 
scheme, accounting for the interference from all other cells. We use tools from stochastic geometry and 
point process theory to analyze the two hop schemes. Besides the results obtained a main contribution 
of the paper is to introduce a mathematical framework that can be used to analyze arbitrary relaying 
schemes. Some of the main contributions of this paper are the analytical techniques introduced for the 
inclusion of the spatial locations of the nodes into the mathematical analysis. 

I. Introduction 

Cellular systems are the most widely deployed wireless systems and provide reliable com- 
munication services to billions around the world. They consist of base stations that serve a 
geographical area called cell. In most of the present cellular systems, the base station (BS) 
communicates directly with the mobile users (MS) in its cell. This single-hop architecture makes 
it is difficult for the BSs to communicate with MSs at the cell boundary because of the distance 
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and the inter-cell interference. So a base station will have to increase its power to maintain the rate 
of transmission. The dead spots problem can be countered by using more base stations, thereby 
increasing the spatial reuse. But increasing the number of base stations can be prohibitively 
expensive or even impossible. The problem can be addressed more effectively by moving away 
from the paradigm of single-hop communication and permitting the base station to communicate 
with mobile stations at the boundary by using the other intermediate MSs in its cell in a sequence 
of hops. Although such multi-hopping requires some significant changes in the present cellular 
system architecture, it may help to effectively combat the dead spots problem, and hence the 
cellular multi-hopping problem is worthy to investigate [1], [2]. In this paper, we analyze the 
benefits of two-hop cellular communication by comparing its performance with a traditional 
single-hop cellular system. A two-hop system, 

• may provide significant benefits over single-hop communication. 

• does not have the implementation complexity of larger number of hops (in terms of routing 
and scheduling). 

When a BS transmits, multiple MSs will be able to receive the information, and hence these 
mobile nodes can help the BS transmit information to the cell edge. Since more than one MS 
can act as a relay, it is not clear how to choose a subset of these relays in a distributed fashion 
so as to reduce the interference and increase the probability of packet delivery. In this paper, we 
analyze simple relay selection schemes and compare their performance with direct transmission. 
We account for the inter-cell interference and the spatial structure of the transmitting nodes in 
the analysis. 

We use methods from stochastic geometry and point process theory to model and study the 
two-hop cellular system. In particular we provide techniques based on probability generating 
functional of a point process to analyze the outage probabilities, and we provide asymptotic 
results for the outage at high SNR and low BS density. The techniques presented in this paper 
can be extended to analyze more complicated relay selection schemes, power control mechanisms 
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and other multi-hop techniques. The major emphasis of the paper is in the methodology and the 
techniques of the analysis rather than the specifics of the communication system. For example 
we concentrate only on two specific relay selection methods although many more methods have 
been proposed in the literature. 

A. Previous work 

The problem of two-hop extensions of cellular system has been studied extensively, and a 
provision for a multi-hop technique has been included in the A-GSM standard [1], [2]. In [3], a 
MS is selected to help the BS depending on the large-scale path-loss on the BS-relay link and 
the relay-destination link. [4] considers a similar problem, but the MSs that can act as relays 
are assumed to be located on a circle around the BS, and the authors provide various power 
allocation schemes and verify their performance by simulations. The present problem is also 
very similar to the problem of opportunistic relay selection. In [5], [6] a detailed analysis of a 
opportunistic two-hop relaying scheme obtaining full diversity order using distributed space-time 
codes has been provided. But a distributed space-time code requires very tight coordination and 
precise signaling among the relays, which increases the overhead and complexity in the system. 
An alternative approach is to choose the best relay, and in opportunistic relaying (OR) [7] a relay 
is chosen so as to maximize the minimum signal-to-noise ratio (SNR) of the source-relay and 
the relay-destination links. In selection cooperation (SC) [8], [9] the relay with maximum relay- 
destination SNR is chosen and has been shown that SC and OR provide a similar diversity order. 
In [5], [7]-[9], distributed relay selection schemes are analyzed and asymptotes of the outage are 
provided for high SNR. The asymptotes provided are functions of the means of fading coefficients 
between the source, relays and the destination. Averaging these results with respect to the spatial 
distribution of the nodes is difficult and hence we use an alternative approach. In our approach 
we model the node locations in a statistical manner and incorporate this information in the 
analysis from the start rather than averaging over the spatial locations at the end. Our emphasis 
is on low-overhead schemes that can readily be implemented. 
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The paper is organized as follows: In Section II the system model is introduced, assumptions 
stated and the metrics used in the paper defined. In Section III the outage probability in the 
direct connection between the BS and its destination is derived. In Sections IV and V the outage 
probability of the two-hop schemes employing different relay selection schemes are analyzed. 
The asymptotic gain of using the two-hop schemes over the direct connection is also studied in 
these sections. In Section VI simulation results are provided and compared to the theory. 

II. System Model 

We assume that the BSs (cell towers) are arranged on a square lattice of density X b . 

The analysis in this paper generalizes in a straightforward manner to any deterministic arrange- 
ment of BS. We assume that n x MSs are available to assist aBSxG $ 6 . More precisely, the 
locations of the mobile stations that assist the base station x form a Poisson point process [10] 
(PPP) of density \ x (y) = rj(y — x). For example choosing rj(y) = l y ([— 1/2, 1/2] 2 ) and 
X b — 1 would lead to a square coverage area for each base station. We use 1 X (A) to denote 
the indicator function of set A. See Figure 1. Observe that it is not necessary for a MS to be 
associated to its nearest BS, i.e., some MSs may be outside the Voronoi cell of their BS. We 
further make the following assumptions: 

1) The average number of MS that each BS serves is finite, i.e., 

N — rj(x)dx < oo. 
Jr 2 

This assumption implies that $ s cannot be homogeneous [10], [11]. 

2) The locations of the mobile users associated with different base stations are independent. 
Since the number of MSs in each cell is Poisson with mean N = E[n x ], each cell is empty 
with probability exp(— N). We shall use /j, to denote the probability that a cell is not empty, i.e., 
pL—1 — exp(— iV). Independent Rayleigh fading is assumed between any pair of nodes and also 
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Fig. 1. Illustration of the cellular system with = 1 and rj(y) — 50 ■ l y ([— 0.25, 0.25] 2 ). So on a average there are 12.5 
MSs per cell. The bold dots represent the BSs and the smaller dots the MSs. The white spaces between the cells may consist of 
other cells which transmit at a different frequency. We may model the case where the neighboring cells use the same frequency 
by choosing rj(y) = l y ([— 0.5, 0.5] 2 ). 



across time, and the power fading coefficient between a node x and node y is denoted by h xy . 
Hence h xy is an exponential random variable with unit mean. The path-loss model is denoted 
by £(x) : M? \ {o} — > M + and is a continuous, positive, non-increasing function of ||x|| that 
satisfies 

£(x)dx < oo, Ve > 0, (1) 



'R 2 \B(o,e) 

where B(a, r) denotes a disc of radius r centered around a. £(x) is usually taken to be a power 
law in one of the forms: 

1) Singular path-loss model: 

2) Non-singular path-loss model: (1 + or min{l, ||x||~ Q }. 

The integrability condition (1) requires a > 2 in all the above models. Assuming simple linear 
receivers and treating interference as noise, the communication between x and y is successful if 

SINR(,, y ,«) = aI + V y KA Xl ' >«■ (2) 
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We also assume 9 > 1 which implies at most one transmitter can connect to a receiver. Here <3> 
is the set of interfering transmitters, p z is the transmission power used by a transmitter located at 
z and a 2 is the the additive white Gaussian noise power at the receiver. We make the following 
assumptions: 

1) In the two-hop schemes that will be analyzed, BSs transmit in the even time slots and the 
MSs transmit in the odd time slots, synchronized across all cells. 

2) Each base station x has an additional mobile station, the destination at r(x) with \\r(x) — 
x\\ = R, to which the BS wants to transmit information. This additional node just receives 
and never transmits. 

3) All the BSs transmit with equal power P. 

Notation: 

• Define 

l(x^y\$) = l(SINR(x, y, $) > 9). 

l(x — > y | <E>) is the indicator random variable that is equal to one if a transmitter at x is 
able to connect to a receiver y when the interfering set is <£>. 

• Define 

$(x) = {?/ G $j : l(x^y\<S> b \ {x})} . 

&(x) is the set of MSs in the cell of BS x to which the BS x is able to connect in the first 
hop (even time slots). 

Metric: Let P d denote the probability that a BS can connect to its destination directly in the 
first hop. Since all BSs are identical 

P d = El(o-r(o) |$ 6 \W)- (3) 

where o denotes the origin (0,0). A BS can connect to multiple MSs in its cell, and these 
connected MS are the potential transmitters in the second hop. In the relay selection methods 
studied in the next section, a subset of these potential transmitters R x C $(x) are selected 
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for each x G $6 to transmit in the next hop. Let the probability that a relay can connect to 
its intended destination (determined by the source to which it connects in the first hop) in the 
second hop be P r , i.e., 

P r = 1 - E J] 1 - l(x - r{o) | \ {x}), (4) 

where ^ = [_} X& ^ B R x is the set of all transmitters in the second hop. Here we are assuming no 
cooperative communication between nodes which have the same information, and hence relays 
belonging to the same cluster also interfere with each other in the second hop. Since 9 > 1, at 
most one transmitter can connect to a receiver and thus 

P r = E]Tl(x^r(o) |M/\{x}), (5) 
and the probability of success for the two-hop scheme is 

P s = l-(1-P d )(l-P r ). 

The BS can potentially transmit in the second hop instead of using the MS as intermediate 
relays. This retransmission scheme will be used as the base reference, and the performance of 
the relay selection schemes will be compared with this retransmission scheme. The gain in using 
the two-hop scheme over the retransmission scheme can be characterized as 

^NM^ (1 _ (i p - ( p -r Prr ^. <« 

where 

P£(R) 



SNR 



a 2 



is the received SNR for the direct transmission. To compare the direct transmission with the 
relay selection scheme, power is allocated across the selected relays in the second hop so that 
the total power is equal to P. Another pertinent metric to capture the performance of the network 
is the diversity gain, defined as 

V bJ SNR^oc log(SNR) 
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From the definition of the diversity and the gain, the following relation follows: 

A(\\ v MG(SNR,A b )) 
d 2 (A 6 ) - d d (X b ) = s hm^ lQg(SNR) , 

where d d is the diversity gain for the single-hop retransmission scheme, and d 2 is the diversity 
gain of the two-hop scheme. From the definition of P s it can be observed that the information 
received in the two time slots is decoded independently. 

In the next sections, we will analyze the success probability P r and the diversity order of the 
relay selection schemes. It is easy to observe that the probability P r of any relay selection scheme 
does not tend to one by increasing the SNR because of the interference caused by transmissions 
in other cells. So to evaluate the asymptotic performance of the system, we scale the BS density 
as 

X b = SHR~ P , (5 > 0. (7) 
As will be evident in the next section, if the signal-to-interference ratio is defined as 

SIR = ^ ^ TTv (S) 



the scaling in (7) translates to 



SIR = 0(SNR^; 



So the system is interference-limited when (5 < 2/a and noise-limited otherwise. Hence the 
scaling in (7) helps us evaluate the performance of the system by varying f3. In practice this 
scaling can be achieved by frequency planning and decreasing the spatial resuse factor. We now 
begin with the analysis of the direct transmission scheme. 
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III. First Hop: Base Station Transmits 

A. Direct Connection 

When the BSs transmit, the inter-cell interference, fading and the noise may cause the trans- 
mission to fail. The probability of direct connection is given by 

P d = El(o -> r(o) | $ fe \ {o}) (9) 

^ >e 



T + E we * fl \{o} h vr(°)^ - r(o)) 

0<7 2 \ 



cxp 



= exp ( -— ) A(r(o)), (10) 



where 



ye* b \{o} 1 + t{x) l ( y ^ 
The following lemma is required to analyze the asymptotics of the success probability. 

Lemma 1: When i(x) = or i(x) = 1/(1 + ||x|| a ), 

l-A(ar) 0C(a) 
lim 77; — = — 

A6-0 A «/2 £( X ) ' 

where 

$(a/2,0)K(a/2 L l/4)-e(a/2,3/4)] 
2 



= • (11) 



£( s > = Sfelo k^-b(k + ~ s i s me generalized Riemann zeta function. 

Proof: We consider the case of i(x) = the other case follows similarly. From the 

definition of A(x) it follows that 

exp^-eeix)- 1 e (y- x )) < A ( x ) < i 1 + ^{x)~ l ^/-z)) -1 - 

ye® b \{o} !/e*b\{°} 

We have 

Y, = £ e (-7T-*) 

y&<3>t\{o} yei?\{o} vv b J 

= K' 2 E *(v-*V>d- 

yeJ?\{o} 
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Dividing both sides by A"/ 2 and taking the limit, the result follows from the definition of the 
Epstein zeta function [12]. ■ 
We have C(3) « 9.03362 and C(4) 6.02681. From the derivation of the above lemma we 
observe that SIR ~ SNR a(3/2 £{R)C(a)- 1 where the definition of SIR is provided in (8). Using 
the above lemma, the asymptotic expansion of for X b = SNR _/3 , (3 ^ 0, at high SNR is 

l-tfSNR" 1 a(3>2 

1-0(1 + C7(a)^(i2)- 1 )SNR- 1 a/5 = 2 (12) 

1 - 6C{a)£{R)- l SNR' a ^ 2 < a/? < 2, 
and the diversity gain of the direct transmission is 

d d (SNR^) =min(l 



2 

So for the direct transmission, (3 < 2/ a corresponds to the interference-limited regime and 
(3 > 2/ a corresponds to the noise-limited regime. From Figure 2 we note that the asymptotes 
in (12) are close to the true 1 — P r even at moderate SNR. In the scaling law provided, observe 
that the distance of the receiver from the BS is fixed. 

B. Properties of the potential relay sets 

In this subsection, the properties of the node set that the BS at the origin is able to connect 
to are analyzed. When the BSs transmit, the interference seen by two MSs is independent. So 
the set of MSs to which the BS at the origin can connect to is an independent thinning of <3> . 
Hence $(o) is also a PPP and since the thinning depends on the position, the resulting process 
is inhomogeneous. Hence the intensity of $(o) is 

8{x) = 7]{x)El{o -> X | $ 6 \ {o}). 

Following a procedure similar to the derivation of (10), the intensity is given by 

6 £(R)\ 



W =^)exp(- — ^JA^). (13) 
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Fig. 2. Outage probability 1 — Pd versus SNR for A;, = SNR~^ with different j3. The system parameters are a = 4, 6 — 1.5, 
r(o) = (0.5, 0.5) and £(x) — (1 + ||x|| 4 ) -1 . The dashed lines are the asymptotes derived in (12). Observe the difference in the 
slopes of the error curve for j3 < 0.5 and j3 > 0.5. 



The average number of MSs which the BS is able to connect to is 

E V l(o->z | $ 6 \{o}) = / 6{x)dx, (14) 

which follows from the Campbell-Mecke theorem [10]. The average distance over which the BS 
at the origin can connect is 

L = ^=g|^ " '' v -—-^ (15) 

E^ t ,l(o-»a!|«i) 

L 2 ||x||(5(a;)dx 

J R2 S{x)dx 

In the second hop, a subset of the MSs which were able to receive information in the first- 
hop transmit. In the next sections we analyze the following two strategies to select a subset 
R x C to transmit in the second hop (odd time slots): 
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• The MS closest to the destination and that has received information in the first hop transmits 
in the second hop. This strategy requires nodes to know their respective locations. 

• The MS with the best channel (fading and path-loss) to the destination that has received 
information in the first hop transmits. This strategy requires the relays to have channel state 
information. 

IV. Method 1 : Nearest Relay to the Destination Transmits 

In this relay selection method, the node x G $(a),a e closest to r(a) is selected to 
transmit in the second hop. To do this each node should know its own location, and each packet 
should have location information about its destination. For a fair comparison with the direct 
transmission scheme, we assume that the selected relay transmits with power SNRcr 2 /£(R). The 
probability of success in this relay selection method is 

Pr = P {Kr(o)i(r) > 6{a 2 + /)) , 

where I is the inter cell interference at r(o), and r is the distance from the relay in the set $(o) 
that is nearest to r(o). More precisely 

J inf x6 $( o) \\x -r(o) ||, |$(o)| >0 

| °°> = o. 

$(o) can be empty because of the following two reasons: 

1) The cell has no MS to begin with. The probability of this happening is 1 — /i. 

2) The BS was not able to connect to any MS in the first time slot. 

For a fair comparison with direct transmission, we condition on the cell at the origin having at 
least one MS to begin with, i.e., n Q > 0. So 

P r | (n a > 0) = P r // _1 . 

Let F a {r, SNR, \ b ) denote the CDF of the first contact distribution of $(o) from r(o). It is given 
by 

F (r,SNR,A 6 ) = l-exp(- f S(x)dx ) . (17) 

V JB(r{p),r) J 
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Observe that F (r, SNR, A b ) is a defective distribution, i.e., F (oo, SNR, X b ) < 1. Let 

/ (r,SNR,A 6 ) = -|-F (r,SNR,A fe ) 



denote the PDF of the first contact distribution. Hence 



f exp (~ 



0£(R) 
SNR£(r) 



E ( exp 



9£(R) 



SNRa 2 £(r)' 

Ti(X b ,r) 



/ (r,SNR,A fe )dr. 



where I is the interference at r(o) caused by transmitters in other cells. Even though J f Q (r, SNR, A b )dr < 
1 the above average is correct since the integrand is zero at r = oo where the remaining mass 
of the first contact distribution lies. We now evaluate Ti(A b , r). Let f a (x), a ^ 0, x E M 2 , denote 
the PDF of the nearest neighbor of r(a) in the set $(a) relative to a, conditioned on the event 
|<&(a)| > 0. We then have 



E 



agZ2\{o} ' 



/a fa) 



4y€(^ + x-r( O ))l(|$(a)|>0) 



dx 



Taking the average with respect to |$(a)| yields 



n 

aSZ 2 \{o} 



E 



l-l(|*(a)|>0) + 



n 1 

aSZ 2 \{o} 



f a (x)l(Ma)\ > 0) 
f a (x)(l - eM~ I S(y)dy)) 



dx 



l + ^(^ + *-r(o))-i 



-dx. 



/ a (x) depends on the geometry of each cell, S(x), r(a), and is easy to calculate once these 
quantities are known. We now calculate the asymptotics of P r and the asymptotic gain. 
Asymptotic gain: In this part we scale the BS density as X b = SNR _/3 . It is easy to observe that 
the average number of MS in each cell that are potential relays, i.e., J S(x)dx, scales as 



J 5{x)dx ~ J r](x)dx - J 



£(x) 



dx — 



6C{a) f 7]{x) 



SNR a/3/2 J £{x) 



-dx. 



(18) 



It can also be verified that 



sup \F a (r, SNR, SNR _/3 ) - F (r, oo, 0)| ^0 
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as SNR — > oo, which implies F Q (r, SNR, SNR ") converges uniformly to F Q (r, oo, 0). Hence we 
can interchange the derivative and the limit in the asymptotic analysis. We have 

/ (r,SNR,SNR-^) = exp (- f 5(x)dx] ^- [ 5(x)dx. 

V JB(r(o),r) J Of J B (r(o),r) 

From (18) and the fact that exp(— x) ~ 1 — x for small x it follows that, 

exp (-/(r)) (&/(r) - «<?(r)) a/3 > 2 
exp (-/(r)) (& /(r) - ^§W(r)) a/5 < 2 



/ (^ SNR, SNR _/3 ) 
where 



= a_ r v(x)_ dx _ r vW dx ! f{r) 

dr J B (r{o),r) Z( X ) JB(r{o),r) H*) dr 

and 

f(r) = / r](x)dx. 

JB(r(o),r) 

The following limit follows similar to the asymptotic analysis of P d 

lim 1 - T i( SNR_/3 ' r ) _ 0C(a)n 
snr^oo SNR _/3a/2 t(r) ' 

where C(a) is given by (11). By some basic algebraic manipulations the asymptotic expansion 
of the error probability with respect to SNR with X b = SNR' 13 , [3 > 0, is 

. pu ^ m J Sf Q °°^P(-f(r)){9(r)+e(r)-^f(r)}dr a(3 > 2 

1 - P2 I [n > 0) ~ ^ " (19) 
1 ^^/o ° ex P(-/W){^)+^W- 1 f/W}dr a/3<2. 

These asymptotes are plotted in Figure 3. From (19) the asymptotic gain is 



lim G(SNR,SNR- /3 ) 

SNR^oo 



^(R)- 1 (j °° exp (-/(r)) {</(r) + €(r)-^/(r)} dr)" 1 a/3 > 2 

^ ^( J R)^(r oo exp(-/(r)){( ? (r)+^(r)- 1 f/(r)}dr)- 1 a/3 < 2. 
Remarks: 

• We observe that the gain is higher in the interference-limited regime than the noise-limited 
regime. This is because in this relay selection method, some of the cells may not be able 
to transmit because they do not contain any MS, which happens with probability 1 — fx. 

• Since the gain does not scale with SNR, the diversity of this scheme is also equal to 
min{l, Pa/2}. See Figure 3 for the error plot obtained by Monte Carlo simulations and the 
above asymptotes obtained theoretically. 
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Fig. 3. Outage probability 1 — P r | (n > 0) versus SNR for At, = SNR~^ for different (3. The system parameters are a = 4, 
8 — 1.5, r(o) = (0.5, 0.5), £(x) — (1 + ||a:|| 4 ) -1 and rj(y) — 51 H ([— 0.5, 0.5] 2 ). The dashed lines are the asymptotes derived in 
(19). The dashed lines are the asymptotes derived in (19) and are approximately equal to 10.351SNR -0 ' 5 (interference-limited) 
and 1.387SNR" 1 (noise-limited). 



V. Method 2: Relay With Best Channel to the Destination 
Transmits(Selection Cooperation). 

In this selection procedure, the fading between a potential relay and the destination is also 
included in the criterion for the relay selection. The relay with the best channel to the destination 
is selected. This method of relay selection is called selection cooperation. In the second hop, 
each relay of the set $(o) can send a channel estimation packet to the destination in an 
orthogonal fashion, and the destination can choose the relay with the best channel. Alternatively, 
if channel reciprocity is assumed, the relays can estimate the channel between themselves and 
the destination when receiving the NACK and use this information to elect the best relay in a 
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distributed fashion. 

As in the previous section we shall find the success probability conditioned on the cell at the 
origin being non-empty, i.e., P r | n Q > 0. As indicated earlier 

P r | (n Q > 0) = fl^Pr. 

Hence we shall first calculate the unconditional probability P r and then multiply it with 
The relay that is selected is mathematically described by 

arg max Ig 4 (o) {h xr{o) £(x -r(o))}. 

The exact analysis of this relay selection in the presence of interference is difficult and hence 
our aim in this section is to obtain the scaling behaviour of G(SNR, SNR _/3 ^ 2 ). Let k denote the 
cardinality of the set $(o). Since the connectivity in the first hop is independent across relays, 
A; is a Poisson random variable with mean 



E[k] 



S(x)dx. 



To make the comparison with the direct transmission easier, we assume that each node transmits 
with power P = SNRcr 2 /£(R). The probability of error is 

1 _ p r = p ( p max {h xr(o) £(x - r(o))} < 6(1 + a 2 ) j , 

where I is the interference at r(o) caused by concurrent transmissions in other cells. Conditioning 
on the point set $(o) we have 



l-P r |$(o) = P P max {h xr{o) £(x - r(o))} < 6(1 + a 2 ) \ $(o) 



E 



If I ex P 

xei>(o) 



6(a 2 + /) 



Ho) 



Pl{x-r(p)), 

Since $(o) is a PPP with intensity function S(x), conditioning on there being k points in the 
set, each node in the set is independently distributed with density k(x) = j /i^dy • Removing 
the conditioning on the locations of $(o), we obtain 



1 - P r | (|$(o)| = k) = E 



-/ 



cxp 



+ 1) 



P£(x - r(o)) 



K,(x)dx 



(20) 
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Using binomial expansion, 

1 - P r | (|$(o)| =k) = E 



E(-D"(*) 

m=0 v 7 



E 



cxp 



fl(a 2 + /) 
'P£(x-r(o)) 







j /t(x)dx 





Hence 1 - P r | (|$(o) | = fc) is equal to 



m=0 

where 



2^( _1 ) \ m J J 2 ^i,---^m)exp ^ jn.(^ 1 ...dx m , 



. . . ,x m ) 



E 



m, 

JJexp 

.6=1 



70 



^ 6 -r(o)) 
= E [exp (-I9q(x™))\ 

= E 



exp -6q(x?) K(a)r(o)((y(a) - r(o))l(\$(a)\ > 0) 

V aeZ 2 / 

where y(a) denotes the location of the selected relay in the cell at a and q(x™) = YlT=i ^( x b ~ 
r(o)) -1 . Let g a (x) denote the PDF of a — x where 

x = arg m&x x ^ (a) {h xr{a) £(x -r(a))}. 

g a (x) is difficult to calculate and is the reason of resorting to asymptotics. Since h y ( ) r ( ) is 
exponential it follows that 

/ \ rr i f g a (y)(l - exp(- J S(x)dx) 

iy(x 1 ,...,X m ) =111-/ m\- lef , a T^W d ^ 



Hence the unconditional probability of error is 



1 — exp I — / 5(x)dx ) 2_] a k 



(/5(x)dx) fe 



fc! 



Asymptotic gain: The above expansion is too unwieldy to yield any asymptotics. We shall use 
(20) to obtain the gain in the high-SNR and low-interference regime. Removing the conditioning 
in (20) we have 

9(a 2 + /) 



1 - P r = Eexp 



-/ 



cxp 



P£(x-r(o)) 



S(x)dx 
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The above result follows from the generating function of a Poisson random variable. Hence 
the required conditional probability is 

9(a 2 + /) 



-/ 



cxp 



P r I (n > 0) = /i" 1 ^1 -Eexp 
An upper bound follows from Jensen's inequality: 

P r | (n > 0) < ( 1 - exp - / Eexp 



Pi{x - r(o)) 



i° 2 + I) 



S(x)dx 



S(x)dx 



P£(x-r(o)) 

Similarly a lower bound can be obtained by using the inequality exp(— x) > 1 — x for the inner 

exp, 



Pr | (n > 0) > fi' 1 1 1 - exp ( - / S(x)dx) Eexp ( / 



-5(x)dx 



P£(x - r(o)) 

To evaluate the upper and lower bounds we observe that we will have to find E[exp(— si)]. By 
a procedure similar to the derivation of v{x\, . . . , x m ): 

Recall that S(x) is equal to 

We now find the asymptotic lower and upper bound when X b = SNR _/3 for large SNR. We first 
observe that 



It is also easy to obtain that 



E[exp(-s/)] ~ 1 - fisC(a). 



After basic algebraic manipulation, it is established that both the upper and the lower bounds 
exhibit the same scaling which is 



Pr | K > 0) 



1 - SNR" 1 [ ±=£ 



SNR-^ 2 (kflocWk 



i + JL 

t(x-r(o)) T t(x) 



r)(x)dx af3 > 2 



U JL 

(x-r(o)) ^ £(x) 



7y(a;)da; a/3 < 2. 



(21) 
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Hence the gain is 



lim G(SNR,SNR" 

SNR^oo 



(x-r(o)) 1 i(x) 



(x-r(o)) ~ £(x) 



ri(x)dx ad > 2 

(22) 

rj(x)dx a(3 < 2 



Hence the diversity of this scheme is 

d 2 (SNR~ /3 ) 



min < 1 



a/3 



In the above analysis we assumed that the cell is non-empty and hence obtained a maximum 
diversity of 1. 



1 -P f | (n o > 0) versus SNR for 9=1 .5, a=4, L=1 using the best relay selection 



P=0 (Monte Carlo) 
P=0.25 (Monte Carlo) 
(3=0.25 Asymptote 
13=0.5 (Monte Carlo) 
13=0.5 Asymptote 
|3=0.75 (Monte Carlo) 
>-p=0.75 Asymptote 




30 40 
SNR dB 



Fig. 4. Outage probability 1 — P r (n > 0) versus SNR for \ b = SNR"' 3 and various f3. The system parameters are a = 4, 
8 — 1.5, z = (0.5,0.5), £(x) — (1 + ||a;|| 4 ) _1 and 77(3/) = 51 y ([— 0.5, 0.5] 2 ). The dashed lines are the asymptotes derived in 
(21) and are approximately equal to 0.812SNR-°' 5 (interference limited) and 0.108SNR -1 (noise limited). 

VI. Simulation Results and Observations 

In this section the gain of the proposed methods over direct transmission is obtained by Monte- 
Carlo simulations. For the purpose of simulation we truncate the BS lattice to X b 1 ^ 2 {—2, —1, 0, 1, 2} 2 , 



November 19, 2009 



DRAFT 



20 



Gain (monte carlo) versus SNR for 8=1 .5, a=4, L=1 




SNRdB 

Fig. 5. G(SNR, SNR _/3 ) versus SNR for various f3. Relay closest to the destination is selected. 

and 9 = 1.5 is used as the decoding threshold. The cells are modeled as squares and the 
destination of each BS is located at a random vertex of the square. The spatial density used is 

V (y) = \ m l y ([-L/2,L/2] 2 ). 

If not specified we use A m = 5 and L — 1. 

()<++> In Figures 3 and 4 the error probability of the schemes employing nearest relay to the 
destination and the best relay are plotted. We observe that the asymptotes obtained from theory 
match perfectly with the simulation results. As predicted by theory, the diversity obtained is 1 
when a/3 > 2 and is equal to af3/2 otherwise. From Figure 5 and 6, it can be seen that the gain 
reaches a constant when SNR — > oo. We observe that the best-relay selection scheme performs 
the best as expected. In Figure 8, we observe that the asymptotic gain increases exponentially 
with \ m because of the (1 — \i)j\x factor in the expression for the asymptotic gain. Setting 
Af, = SNR _/3 reduces the spatial reuse factor as the SNR increases. The effective throughput 
density of the network is equal to P 2 log(l + f9)SNR~ /3 and the maximum of this throughput 
density is the transmission capacity [13]. In Figure 9, we plot (P 2 | n > 0)log(l + #)SNR -/3 
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Fig. 6. G(SNR, SNR ) versus SNR for various (3. Relay with the best channel to the destination is selected. 



Gain versus X m , r)(y)=X. m 1 y ([-L/2 L/2] ), 8=1.5 for the nearest relay selection 



4.5 
4 
3.5 




- |> » t 




Interference limited L=1 
Noise limited L=1 
Interference limited L=2 
Noise limited L=2 



0.5 L 




5 10 
Mobile Intensity X 



15 



Fig. 7. Asymptotic gain versus A m where A m is the intensity in rj(y) — \ m l y ([— L/2, L/2] 2 ), t{x) = a , 
z = (-L/2,L/2). 



= 1.5 and 
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Gain versus A. , t](y)=X 1 ([-L/2 L/2] 2 ), 8=1 .5 for the best relay selection 

io 3 m 




Fig. 8. Asymptotic gain versus A m where A m is the intensity in r)(y) = A m l !/ ([— i/2, L/2] 2 ), i{x) = a , 9 = 1.5 and 
z = (-L/2,L/2). 

versus SNR for various (3. We observe that for each SNR there is a (3 that maximizes the 
throughput density, and that as SNR — > oo, the maximizing f3 tends to 0, which is intuitive. The 
figure indicates that a throughput density of ss 0.1bps/m 2 is achieved at low SNR, and that it 
increases with SNR. 

VII. Conclusions 

In this paper we have analyzed the outage in a two-hop cellular system under consideration 
of all the node location statistics. Outage results were provided for two relay selection schemes, 
namely nearest-relay selection and best-relay selection. We observed that the diversity obtained 
is min{l, a{3/2} where a is the path-loss exponent, when the density of the base stations scale 
as X b = SNR _/3 (alternatively SIR = 0(SNR a/3/2 )). From this result we can infer that the system 
is noise-limited (even for high SNR) when a(3 > 2 and interference-limited otherwise. The 
asymptotic outage gain of the two-hop system over direct transmission takes only two values as 
a function of (3 depending on the relay selection scheme. The gain in selecting a relay with the 
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Fig. 9. (P'2 \ n a > 0) log 2 (l + #)SNR * versus SNR for various (3. The best relay selection scheme is used. 



best channel over a direct transmission increases exponentially with the density of the available 
relays. The gain also increases with increasing source-destination distance. From simulations 
we conclude that the gain in selecting the best relay outweighs the overhead in estimating the 
fading coefficients between the relays and the destination as compared to the near-relay selection 
method. The techniques introduced in this paper can be extended for the spatial analysis of other 
relay selection schemes. 
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